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SECTION  1 


INTRODUCTION 

Several  American  and  Japanese  institutions  are  engaged  in  research 
efforts  to  produce  and  validate  a  revised  dosimetry  system  for  radiobiological 
analysis  of  the  survivors  of  the  Hiroshima  and  Nagasaki  A-bombs.  *  The 
revised  system  will  be  produced  by  modern  calculational  techniques  as  dis¬ 
tinguished  from  the  dosimetry  system  presently  used  (called  T65D)  which 
relied  mainly  on  Nevada  Nuclear  Weapons  Tests  and  the  BREN  (Bare  Reactor 
Experiment  f^evada)  experiments.  In  situ  data,  the  Nevada  Weapon  tests,  the 
BREN  experiments,  and  other  appropriate  experimental  data  are  being  used  to 
validate  the  techniques  and  computer  codes  used  for  calculations  of  the 
weapon-produced  radiation  sources,  the  air-over-ground  radiation  transport, 
the  house  shielding,  and  the  in-body  radiation  transport.^ 

Work  on  calculating  the  free-in-air  (FIA)  radiation  environment  is 
well  advanced  with  detailed  two-dimensional  calculations  of  the  radiation  out¬ 
put  of  the  asymnetric  LITTLE  BOY  device  exploded  over  Hiroshima  completed  by 
2 

Whalen  and  with  subsequent  calculations  of  the  prompt  radiation  environment 

2 

by  Pace  and  the  delayed  radiation  by  Scott.  An  independent  calculation  of 
the  FIA  environment  at  both  Hiroshima  and  Nagasaki  from  older,  one-dimen- 
sional  weapon  output  radiation  computations  by  Preeg  for  the  A-bombs  has 
been  reported  by  Lcewe.^ 

The  FIA  radiation  environment  can  be  significantly  perturbed  by  the 
house  occupied  by  the  survivor  and  by  nearby  neighboring  houses.  Compared  to 
the  FIA  environment,  the  exposure  environment  at  the  survivor  location  may  be 
reduced  by  factors  from  two  to  four,  be  shifted  in  energy  and  have  a  reduced 
neutron  to  gamma-ray  ratio.  We  are  engaged  in  a  program  to  calculate  the 


house  shielding  effects  on  the  exposure  radiation  environment  and  to  in¬ 
corporate  these  effects  into  the  revised  dosimetry  system. 


As  a  first  step  in  this  research  effort,  we  have  validated  our 
method  by  calculating  the  BREN  house  transmission  experiments®  which  were 
used  to  derive  the  house  shielding  for  the  T65D  system.  Although  we  believe 
that  the  T65D  house  radiation  transmission  factors  are  in  error  and  need 
revision,^  the  BREN  experiments  can  be  used  to  validate  the  calculational  pro 
cedures.  The  same  techniques  can  be  used  to  determine  transmission  factors 
with  more  realistic  models  of  the  A-bomb  events  than  could  be  achieved  in 
experiments  such  as  BREN.  This  paper  presents  our  work  to  compare  calculated 
house  transmission  factors  with  reported  measured  data  from  the  BREN  experi¬ 
ments  . 
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SECTION  2 
BACKGROUND 

Detailed  descriptions  of  the  location  and  position  at  the  time  of 
bomb  detonation  (ATB)  of  the  A-bomb  survivors  in  the  Radiation  Effects  Re- 

g 

search  Foundation  (RERF)  data  bases,  such  as  Life  Span  Study  (LSS),  have 
been  compiled  from  interviews  and  placed  in  shielding  files.  Nine  categories 
were  used  to  describe  the  shielding  situation  of  the  survivor.  The  number  of 
subjects  in  each  category  for  the  proximal ly  exposed  group  (less  than  1600-m 
ground  range  at  Hiroshima  and  less  than  2000-m  at  Nagasaki)  is  given  in  Table 
1.  Over  half  of  the  survivors  in  the  data  base  were  inside  residential 
houses  ATB. 

In  the  T65D  system,  the  perturbation  of  the  FIA  radiation  by  res¬ 
idential  houses  was  estimated  by  transmission  factors  for  neutrons  and  gamma 
rays  which  multiplied  the  FIA  kermas  to  give  the  exposure  neutron  and  gamma- 

3 

ray  kerma .  Later,  organ  dose  factors  were  developed  to  account  for  body 

9 

shielding.  The  transmission  factors  for  individual  survivors  were  calculated 
by  a  parametric  model  called  the  "Nine  Parameter  Formula"  (NPF)  based  on  a 
regression  analysis  of  the  BREN  house  transmission  experiments.^  The  nine 
parameters  were  derived  for  individual  survivors  from  examination  of  his/her 
shielding  file  and  coded  for  computer  access.^® 

The  BREN  house  transmission  experiments  consisted  of  exposing  ra¬ 
diation  analogs  of  Japanese  residential  houses  to  radiation  from  the  bare 
Health  Physics  Research  Reactor  (HPRR)  and  a  ^*^Co  source.®  The  sources 
were  suspended  up  to  1125  feet  from  a  tower  and  the  houses  were  placed  at 
ground  ranges  from  750  to  1200  yards  from  the  tower.®  The  transmission  factors 
reported  were  the  ratio  of  the  neutron  and  gamma-rays  kermas  measured  inside 

®In  the  following,  we  will  use  the  dimensional  specifications  (inches,  feet,  yards) 
aooearing  in  the  BREN  experiment  reports  for  clarity  and  ease  of  comparison. 
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the  houses  at  various  locations  to  kemas  measured  at  the  same  ground  range  at 
locations  free  from  local  obstructions.  The  neutron  measurements  were  made 
with  a  recoil  proton  proportional  counter  and  the  gamma-ray  measurements  were 
performed  with  Geiger-Muel ler  tubes. 

With  this  experimental  procedure,  the  ganma-ray  kerma  measured  in¬ 
side  the  house  includes  garrma  rays  produced  by  neutron  inelastic  and  capture 
reactions  with  the  house  materials.  Thus,  the  gamma-ray  transmission  factors 
are  dependent  on  the  local  neutron  environment.  For  the  nominal  configura¬ 
tion  at  BREN  (HPRR  at  1125  ft  height,  houses  at  750  yds)  the  calculated 
neutron  to  ganma-ray  kerma  ratio  is  nearly  three,  while  at  Hiroshima  and 
Nagasaki  the  same  ratio  is  less  than  0.1  and  0.01,  respectively,  at  ground 
ranges  where  survivors  in  the  data  bases  were  present  (beyond  about  900  m). 
This  problem  was  first  examined  by  Marcum  who  concluded  that  the  gamma-ray 
transmission  factors  used  in  T650  based  on  BREN  data  were  too  high  by  a 
factor  of  1.6^’^^  because  of  the  higher  house  gamma  ray  contribution. 

Although  the  BREN  experiments  did  not  accurately  replicate  the 
radiation  environments  at  Hiroshima  and  Nagasaki,  they  are  valuable  experi¬ 
ments  for  validating  calculational  procedures  which  can  then  be  applied  to 
the  Japanese  house  shielding  problem.  Our  method  to  calculate  the  house 
shielding  uses  the  Monte  Carlo  technique  in  the  adjoint  mode  to  compute  the 
energy-angular  radiation  response  of  the  house  cluster  surrounding  the 
detector  location.  This  response  function  can  then  be  coupled  to  the  FIA 
energy-angular  environment  to  estimate  the  exposure  environment. 
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SECTION  3 


THEORY 

The  radiation  environment  inside  a  small  structure  far  removed  from 
a  source  of  radiation  can  be  determined  by  partitioning  the  problem  into  sep¬ 
arate  transport  calculations  that  are  then  coupled  to  provide  the  overall 
solution.  This  partition  consists  of  a  calculation  of  the  global  radiation 
environment  in  simplified  geometry  coupled  to  an  appropriate  radiation  trans¬ 
port  calculation  in  a  detailed  geometry  model  of  the  structure.  The  first 
calculation  is  a  discrete  ordinates  numerical  solution  to  the  Boltzmann 
transport  equation  (BTE)  in  two-dimensional  air  over  flat  ground  geometry; 
the  latter  calculation  is  often  a  Monte  Carlo  solution  to  the  adjoint  BTE  in 
a  three-dimensional  model  of  the  structure. 

The  theoretical  basis  for  this  solution  technique  is  the  fact  that 
the  radiation  environment  inside  a  volume  enclosed  by  a  non-reentrant  surface 
can  be  completely  characterized  by  the  knowledge  of  the  incoming  particle  flux 
at  the  surface  and  the  radiation  sources  within  the  volume.  Since  the  struc¬ 
tures  here  have  no  internal  sources  of  radiation,  we  can  determine  the  inte¬ 
rior  radiation  environment  if  we  know  the  incoming  particle  flux  on  a  non-re¬ 
entrant  surface,  e.g.,  sphere  or  cylinder  enclosing  the  structure. 

The  theory  for  this  procedure  has  been  rigorously  developed  and 
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adequately  explained  by  many  authors,  notably  Hoffmann,  et  al.  We  provide 
here  only  a  summary  to  maintain  continuity  for  the  reader.  Following  Hoffman, 
the  rigorous  solution  to  the  complete  problem  (including  both  source  and 
structure  geometry)  is  given  by 


H4>(F)  =  S(T^) 

6 


(1) 


where  H  is  the  Boltzmann  operator;  4>(P)  is  the  flux  in  six-dimensional  phase 
space  denoted  by  P”;  and  S(F)  is  the  source.  The  result  of  interest  X  (i.e., 
some  observable  quantity)  involves  integration  of  a  known  flux  response  func 
tion  R(F)  with  the  particle  flux. 


X  =  yi(F)R(P)dP  . 


Now,  we  assume  that  the  source  S(F)  is  zero  in  the  region  of  the 
structure  and  the  response  function  R(F)  is  zero  outside  the  region  of  the 
structure.  Let  r^  denote  the  vector  of  the  locus  of  points  describing  a 
nonreentrant  surface  enclosing  volume  V  containing  the  structure.  Let  F^ 
denote  phase  space  elements  in  7.  We  require 


S(P,)  =  0 


R(F)=0.  F/F^  . 


In  this  case,  the  rigorous  solution  for  the  flux  in  V  can  be  written  as 


H'})(Fy)  =  -(n*n)(J){Fy)6(r-rg)  ;  n*fl  <  0,  (5) 


where  7  is  the  outward  directed  normal  to  the  non-reentrant  surface,  and  the 
spatial  and  angular  elements  of  F  are  denoted  by  F  and  ?i,  respectively.  We 
have  replaced  Eq.  (1)  for  the  flux  in  7  by  an  equivalent  expression  with 
Eq.  (5). 


The  well-known  property  of  the  Boltzmann  operator 


( 


J  <J>*(P^)H4)(P^)dP^  =  f  (t>(P^)H*4>*(Py)dPy 


+  y*  (n*fi)^(P^)<j)*(Fy)6(r-rj)dPy  (6) 


reduces  to 


^  =  /<KPv)R(Rv)dP, 


./* 


ft)4>(Py)4>*(?^y)6(r-r5)dPy 


by  noting  that  the  LHS  of  Eq.  (6)  is  zero  since  H(})(F^)  =  S(Py)  =  0  and  by 


substituting 


H*0*(Py)  =  R(Py)  . 


Here,  H*  is  the  adjoint  Boltzmann  operator.  By  virtue  of  Eq.  (5),  the  result. 


A,  cannot  depend  on  the  outward  directed  flux  on  the  surface  of  v.  This 


suggests  the  boundary  condition  for  the  solution  for  <))*{Py)  as 


(n*fi)^>*{P  )6(r-r  )  =  0,  n»fi  >  0 

V  « 


Thus,  by  solving  Eq.  (8)  for  with  boundary  condition  given  by  Eq.  (9) 


(vacuum  boundary),  the  result.  A,  can  be  obtained  via 


vv 


ol 


J\ 


Our  solution,  at  this  point,  has  no  approximations.  But  the  flux 
in  Eq.  (10)  contains  the  structure  in  the  problem.  We  now  solve  for 


a  new  flux 


Va^^^  =  S(P) 


for  the  same  source  but  in  a  simpler,  approximate  geometry  without  the  details 
of  the  structure  present,  as  embodied  in  the  operator  We  then  form  the 

fl 

approximate  result 


•fi)<t'a(^v)<l>*(^„)6{f-rg)dPy,  n*0  <  0, 


with  an  error 


■  7‘" 


Note  that  4i*(Py)  is  still  the  rigorous  solution  in  P^. 

To  summarize,  we  calculate  the  effect  of  interest  inside  a  complex 
structure  in  air-over-ground  geometry  by  performing  a  surface  integral  of  the 
adjoint  flux  computed  in  the  complex  geometry  with  the  incoming  surface  flux 
from  the  source  computed  in  a  simplified  geometry.  The  error  of  this  approxi¬ 
mate  approach  is  directly  associated  with  the  difference  in  the  incoming  flux 
at  the  surface  between  the  simplified  geometry  and  the  complex  geometry. 


Thus,  assuming  an  accurate  adjoint  solution  <>i*(Py)  can  be  obtained,  assuming 
the  incoming  flux  on  the  surface  surrounding  the  structure  is  not  significantly 
perturbed  by  its  presence,  and  assuming  the  calculation  of  the  flux, 
in  the  simplified  geometry  model  is  accurate;  the  solution  given -by  Eq.  (12) 
should  be  a  reasonable  •jstimate  of  the  desired  result  in  Eq.  (2). 
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SECTION  4 

APPLICATION  OF  THE  THEORY 

A  code  system  called  VCS  (Vehicle  £ode  System)^^  was  developed  to 
calculate  6(7-7^)  and  perform  the  numerical  intergration  of  .Eq.  (12) 

with  the  forward  fluxes,  from  a  calculation  with  the  two-dimensional 

d  V 
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discrete  ordinates  code,  DOT.  This  code  system  was  originally  designed  to 
calculate  radiation  protection  factors  for  military  vehicles  such  as  tanks 
and  armored  personnel  carriers.  We  have  adapted  the  VCS  code  with  some  major 
and  minor  modifications  for  the  house  shielding  analyses.  In  the  following, 
we  discuss  how  our  modified  VCS  code  system  is  used  to  estimate  the  radiation 
environments  in  structures. 

The  calculation  of  the  forward  flux,  4'-(7,),  is  performed  with  the 
DOT  code.  The  DOT  code  uses  multigroup,  multitable  cross  sections  coupling 
neutron  transport,  gamma-ray  production  from  inelastic  and  capture  reactions, 
and  gamna-ray  transport.  The  method  utilizes  discrete  ordinates  in  two- 
dimensional  geometries  (spherical,  cylindrical,  and  slab).  The  cylindrical 
geometry  is  used  for  VCS  calculations  consisting  of  air  over  a  flat  ground 
plane  with  the  source  located  above  the  ground  on  the  axis  of  symmetry.  This 
simplified  two  material  geometry  is  adequate  for  calculations  of  the  BREN 
reactor  and  ®^Co  source  suspended  above  the  desert  in  Nevada;  however,  cal¬ 
culations  for  built-up  areas  such  as  Hiroshima  and  Nagasaki  may  require  a 
mush  model  of  buildings  to  account  for  the  distal  effects  of  the  built-up 
urban  areas  on  localized  radiation  environments. 

The  DOT  code  for  cylindrical  geometry  provides  a  computer  file  of 


the  angular  flux  quantity 


<1)3(2, r,n,<J^,E) 


where  the  spatial  coordinates  (z,r)  and  angular  coordinates  (n,)!))  are  defined 
in  Figure  1.  Note  that  n  refers  to  the  cosine  of  angle  of  the  angular  direc¬ 
tion  with  the  z  axis  and  ii  is  the  angle  of  the  projection  in  the  x-y  plane 
with  respect  to  the  radial  vector.  (The  subscript  "a"  for  approximate  is 
retained  to  note  that  the  calculation  does  not  include  the  specific  structure 
in  which  the  radiation  environment  is  desired.)  The  flux  is  provided  for  dis¬ 
crete  mesh  intervals  in  z,r  space,  in  11,4)  angular  space,  and  for  energy  groups 
in  E.  Let  co.  .  be  the  fraction  of  4tt  angular  space  ascribed  to  the  discrete 
ordinate  (r].,<p..).  The  scalar  flux  can  be  calculated  as 

I  I J 


43(2,r,E)  =  I  I 

i  j 


(14) 


The  units  on  the  angular  flux  from  DOT  are  particles  per  cm'^  per  fraction  of 
47r  space.  The  VISA  code  of  VCS  is  used  to  access  the  DOT  angular  flux  file 
and  retrieve  only  flux  values  near  the  air/ground  interface  required  for  the 
coupling  calculation  of  Eq.  (12). 

The  adjoint  quantity,  <i)*(F^)6(r-7j) ,  is  computed  with  an  adaptation 
of  the  MORSE  Monte  Carlo  code^^  called  MORSE/VCS.  The  MORSE  code  uses  multi¬ 
group,  multi  table  cross  sections  and  operates  in  both  the  forward  and  adjoint 
modes.  Special  user  routines  can  be  written  and  interfaced  with  MORSE  to  handle 
source  specifications  and  flux  estimation.  Transport  calculations  can  be  car¬ 
ried  out  in  complex  geometrical  configurations  modeled  with  the  combinatorial 
geometry  package  in  MORSE. 


The  primary  modification  to  MORSE  in  the  MORSE/VCS  code  is  the  in¬ 
clusion  in  the  code  of  routines  to  write  a  particle  leakage  file  containing 
the  source  energy  group,  the  leakage  energy  group,  the  leakage  direction 
cosines,  and  the  "z"  coordinate  (height  above  the  ground  plane).  A 
particle  "leaks"  in  the  adjoint  MORSE/VCS  calculation  when  it  crosses  the 
vacuum  boundary  surface  that  defines  the  volume  including  the  structure,  in 
accordance  with  the  boundary  condition  specified  in  Eq.  (9).  This  volume  is 
a  right  circular  cylinder  with  the  axis  normal  to  the  ground.  The  cylinder 
entirely  encloses  the  object  of  interest  and  extends  a  few  tens  of  centi¬ 
meters  into  the  ground. 

Minor  modifications  that  have  been  made  for  these  house  shielding 
calculations  include  the  provision  to  write  on  the  leakage  file  the  entire 
source  and  leakage  phase  space  parameters.  This  will  enable  eventual  coupl¬ 
ing  to  man  phantom  calculations  to  estimate  body  shielding.  The  geometrical 
surface  (e.g.,  roof,  right  side  window,  etc.)  that  the  particle  crossed  when 
escaping  into  the  surrounding  air  or  ground  is  also  written  for  later  inter¬ 
pretation  of  the  results. 

The  major  modification  to  the  code  was  the  inclusion  of  adjoint  in¬ 
group  biasing  developed  by  Scott^^  to  improve  the  statistical  precision  of 
adjoint  coupled  calculations  where  gamma  rays  become  thermal  neutrons,  and 
then  the  thermal  neutrons  upscatter  to  higher  energies.  If  these  calculations 
are  performed  in  the  standard  manner,  large  variations  in  particle  weights 
occur  (hence,  large  standard  deviations  in  results)  becuase  the  adjoint  ther¬ 
mal  neutrons  are  "forced"  to  upscatter  with  a  large  fair  game  weight  correction. 


The  adjoint  calculation  is  performed  with  a  detailed  mathematical 
model  of  the  structure  using  the  combinatorial  geometry  package  available 
with  MORSE.  The  model  contains  the  ground  and  air  surrounding  the  struc¬ 
ture  inside  the  leakage  surfaces  of  the  cylinder.  The  adjoint  sfmulation 
commences  at  the  detector  location  in  the  structure.  As  the  adjoint  particle 
interacts  with  the  materials  in  the  building,  air  and/or  ground,  it  gains 
energy  and  gamma  rays  produce  adjoint  secondary  neutrons.  The  selection  of 
scattered  energy  and  angle  uses  the  same  procedures  and  routines  as  the  forward 
calculational  procedure,  except  the  reaction  cross  sections  for  adjoint  cal¬ 
culations  are  inverted  from  the  forward  set.  When  a  particle  leaks  from  the 
system,  the  relevant  parameters  are  written  to  the  leakage  file  and  a  new 
history  is  started. 

When  the  calculation  is  completed,  this  leakage  file  contains  the 
individual  statistical  weights  of  all  the  leaking  particles  as  a  function  of 
source  and  leakage  variables.  This  may  be  written  as 

W  t  Rj  •  fig  .^5 ) 

where  W  is  the  statistical  weight  of  the  particle,  (R|_,fi^_,E|^)  and 

are  leakage  and  source  position,  angle,  and  energy  phase  space  parameters, 

respectively.  In  the  following  we  also  employ  the  notation 

to  decompose  the  spatial  vector  into  individual  orthogonal  components. 


The  DRC  code  of  VCS  reads  the  DOT  forward  flux  tape  prepared  by 
VISA  and  the  adjoint  leakage  file  provided  by  MORSE/VCS,  and  then  performs 


the  numerical  integration  of  Eq.  (12). 

Because  the  structure,  in  principle,  can  be  oriented  in  any  direc¬ 
tion  with  respect  to  the  source,  the  angular  coordinates  of  each  leakage 
particle  are  transformed  to  the  desired  orientation.  This  transformation 
with  a  matrix  derived  from  Eulerian  angle  rotation  is  given  by 


=  T  •  (16) 

=  T  .  \  (17) 

where  T  is  the  transformation  matrix.  The  standard  VCS  code  only  permits  rota¬ 
tion  of  the  structure  about  an  axis  parallel  to  the  Z  direction.  We  have  added 
the  generalized  rotation  to  examine  effects  of  small  "tilts"  in  the  structure, 
even  though  the  ground  planes  of  the  DOT  and  adjoint  calculations  would  not 
precisely  align. 

With  the  transformed  angle  coordinates,  the  angular  bin  of  the  DOT 
quadrature  containing  the  leakage  direction  is  found.  First,  the  "eta"  level 
is  determined  based  on  the  value 

n  =  z  •  ,  (18) 

where  z  is  the  unit  vector  normal  to  the  ground  plane.  The  relative  weight 
(or  fraction  of  4tt  space)  associated  with  each  polar  eta  level  is  calculated 
by  DRC  via 
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and  "i"  index  (eta  level)  is  found  such  that 


i-1  i 

[  w. ,  <  i  (n  +  1)  <  I  co.i  . 
i'=l  '  i-=l  ' 


The  phi  level  ("j"  index)  is  based  on  the  value: 


cos  =  r  •  •  (1  -  ) 


(20) 


(21) 


(where  r  is  the  unit  vector  in  the  radial  direction)  and  is  found  via 


J'l  ,  ,  J 

I  W.  .  ,  <  y  (cos  tl'.+l)  <  ~  I  U). 
j'=l  2  1  IJ 


(22) 


This  procedure  determines  the  discrete  ordinate  .  =  (n.,  4'..)  for  the  DOT 

I J  1  1 J 

angular  flux. 

The  estimate  of  exposure  flux  by  Eq.  (12)  is  then  given  by 


^(Ts.n’s.Es)  =  <  <f»3(zL.rj,P.j.EL)-W(RL,n’L.EL,  Rs.^^'s.E5)>  (23) 

where  r^  the  ground  range  of  the  in-structure  detector  from  the  source  in  the 
DOT  calculation  and  where  the  brackets  <>  indicate  the  average  over  the  leak¬ 
age  variables  for  all  histories.  Note  that  this  estimate  only  depends  on  the 
"z"  spatial  coordinate  of  the  adjoint  leakage.  The  forward  DOT  flux  for  the 
coupling  is  only  defined  at  the  radius  of  the  detector  location  r^.  This 
assumes  that  the  variation  in  flux  with  radius  is  small  compared  to  the  radial 
extent  of  the  structure. 

17 


By  preserving  the  full  phase  space  parameters  of  the  adjoint  source. 


the  full  energy  angular  distribution  of  the  exposure  flux  can  be  calculated. 
This  can  be  used  for  detailed  calculations  of  radiation  transport  in  man 
phantoms  using  the  same  surface  integral  technique  (Eq.  12)  with -an  addi¬ 
tional  adjoint  calculation  in  the  man  phantom  starting  at  the  organ  location. 
This  facility  to  "nest"  calculations  using  successive  adjoint  calculations 
was  incorporated  by  us  into  the  standard  version  of  VCS,  This  is  done  by 
defining  the  angular  flux  in  the  same  angle  bins  as  the  DOT  quadrature,  so 
that  the  flux  calculated  from  a  coupling  procedure  is  in  the  same  format  as 
the  VISA  processed  DOT  angular  flux. 

The  exposure  flux  inside  the  structure  can  be  used  to  calculate 
detector  responses  and  kermas  via 


X(r3) 


/• 


=  /<^(r5,ii'5,E5)R(Ej)dO'5dE5 


(24) 


where  R(E^)  is  the  energy-dependent  response.  Note  that  if  the  energy  angular 
dependent  response  is  known,  it  can  be  integrated  with  the  angular  dependence 
of  the  exposure  flux. 


In  addition  to  the  approximation  introduced  by  the  theory,  the 
practical  application  with  the  VCS  methodology  introduces  approximations  and 
sources  of  uncertainty.  There  exist  the  usual  sources  of  uncertainty  asso¬ 
ciated  with  any  calculational  technique,  including  the  knowledge  of  the  phys¬ 
ical  variables  associated  with  configuration  (e.g.,  the  air  density,  construc¬ 
tion  materials,  structure  layout);  the  basic  data  input  to  the  calculation, 

’  such  as  nuclear  cross  sections;  and  modeling  of  the  configuration  and  the 

j 
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physics  data  (e.g.,  number  of  geometry  regions,  energy  group  structure, 
quadrature  order). 

There  are  also  sources  of  uncertainty  peculiar  to  the  surface 
integral  coupling  method.  One  of  the  most  important  is  the  coupling  to  DOT 
fluxes  at  a  single  radial  point  corresponding  to  the  ground  range  of  the 
detector.  Although  VCS  performs  the  surface  integral  as  a  function  of  a  z 
coordinate,  the  code  assumes  the  flux  does  not  vary  in  the  radial  direction 
across  the  leakage  surface.  To  analyze  this  assumption,  we  write  the  radial 
variation  of  the  flux  approximately  as 

(25) 

r 

for  some  appropriate  relaxation  length,  L.  This  is  a  good  representation 
overa  limited  rangeat  least  and  was  used,  for  example,  in  the  T65D  system  to 
provide  the  doses  at  all  ground  ranges  used  by  RERF.  For  small  Lr,  we  can 
write 


<t)(r)  -  <|)(r+Ar)  „  Z_\ 

m  *  r  j 


(26) 


The  relaxation  length  is  on  the  order  of  300  m  for  many  sources  in  air-over¬ 
ground  and  choosing  r  to  be  the  same,  we  estimate  as  a  near  worst  case  (since 
r  is  usually  larger),  that  the  flux  changes  by  10%  in  10  m.  This  indicates 
that  the  assumption  is  reasonable  for  small  structures;  however,  the  method 
may  have  to  be  modified  for  problems  involving  a  cluster  of  houses  for  which 
there  may  be  unacceptable  variation  in  the  flux  across  the  larger  volume 
enclosed  by  the  leakage  surface. 
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Another  issue  more  difficult  to  analyze  is  the  adequacy  of  the 
resolution  of  the  angular  flux  in  the  DOT  quadrature.  The  VCS  procedure 
assumes  the  flux  is  constant  within  the  solid  angle  defined  by  the  quadrature 
boundaries.  Furthermore,  the  uncollided  flux  is  "smeared"  over  the  most  for¬ 
ward  directed  quadrature  interval.  The  concern  may  be  that  streaming  paths, 
such  as  windows,  may  receive  too  much  or  too  little  incident  flux,  depending 
on  its  relation  to  the  quadrature  boundaries  and  the  "true"  angular  variation 
of  the  flux. 
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SECTION  5 


COMPARISON  WITH  THE  BREN  EXPERIMENTS 

5.1  BREN  HOUSE  TRANSMISSION  MEASUREMENTS. 

The  Operation  BREN  experiments  took  place  at  the  Nevada  -Test  Site 
(NTS)  in  the  spring  of  1962.^^  The  experiments  consisted  of  suspending  a 
small  bare  reactor  and  a  ®^Co  source  from  a  tower  at  heights  of  up  to  1125 
feet  and  making  measurements  of  free  field  and  house  perturbed  neutron 
and  gamma-ray  kermas.  The  ratio  of  the  kermas  measured  at  detectors  in  (and 
around)  the  houses  to  free  field  kerma  was  reported  as  the  house  transmission 
factor  for  the  detector  location. 

The  Health  Physics  Research  Reactor  (HPRR)  was  used  for  the  neutron 
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(and  secondary  gamma  ray)  source.  The  reactor  consisted  of  a  highly  enriched 
uranium-molybdenum  alloy  metal  annulus  with  an  outside  diameter  of  0.203  m 
and  inside  diameter  of  0.05  m.  Details  of  the  reactor  configuration  are 
found  in  Ref.  18.  The  small  size  of  the  reactor  and  lack  of  reflector  pro¬ 
vides  a  hard,  nearly  unmoderated  fission  spectrum  of  neutrons.  The  ^®Co 
source  used  at  BREN  contained  1200  Ci  at  time  of  manufacture  which  had 
decayed  to  800  Ci  at  the  time  of  the  experiments. 

The  free  field  ganna-ray  kermas  produced  by  the  HPRR  and  ®*^Co  source 
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were  measured  with  the  "Phil"  detector  described  by  Wagner  and  Hurst.  It 
consists  of  a  small  Geiger-Mueller  tube  which  was  enclosed  by  tin  and  lead 
sheaths  to  provide  (according  to  the  developers)  a  uniform  kerma  per  count 
as  a  function  of  incident  gamna-ray  energy  from  100  keV  to  10  MeV.  However, 
at  the  time  of  the  BREN  experiments,  the  highest  reported  calibration  energy 
was  the  ^®Co  gamma-ray  lines  (1.17,  1.33  MeV).  When  a  neutron  field  was  pre¬ 
sent,  the  Phil  was  further  enclosed  by  a  ®Li  shield  to  reduce  the  neutron 


response  from  thermal  neutron  capture.  It  was  claimed  that  the  Phil  had  no 
fast  neutron  response.  To  improve  counting  statistics,  measurements  inside 
the  houses  were  performed  with  a  larger  Geiger-Mueller  detector.  Details  of 
the  calibration  or  encapsulation  materials  for  this  detector,  however,  were 

never  reported  in  the  open  literature.  The  neutron  kerma  measurements  were 
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made  with  the  RADSAN  proportional  counter. 


The  houses  constructed  at  NTS  for  the  BREN  experiments  were  designed 

to  be  radiation  analogs  of  typical  Japanese  residential  houses  at  Hiroshima 

and  Nagasaki  ATB.^  Three  models  were  built  designated  as  House  A  (a  typical 

one-story  family  house);  House  B  (a  larger  two-story  residence),  and  House  C 

(a  small  one-story  tenement  house).  The  houses  were  constructed  full  scale 

to  the  very  uniform  dimensions  of  Japanese  dwellings.  In  place  of  Japanese 

construction  materials,  however,  the  designers  used  a  cement  asbestos  board 

produced  by  the  Johns  Mansville  Corporation  under  the  trade  name  TRANSITE, 

which  was  used  in  many  houses  in  East  Tennessee  near  the  Oak  Ridge  National 

Laboratory  (ORNL).  The  use  of  cement  asbestos  board  was  justified  by  the 

results  of  direct  beam  radiation  measurements  at  ORNL  which  showed  the  ratio 

of  PuBe  neutron  attenuation  to  ®®Co  attenuation  similar  to  the  ratio  for  a 
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Japanese  house  wall  section.  The  exterior  walls  consisted  of  1-inch  and 
3/4-inch  thick  sheets  butted  together  and  placed  in  slots  on  wood-framing 
members.  The  interior  walls  were  actually  thicker  and  consisted  of  two 
1-inch  sheets.  The  roof  had  the  same  total  thickness  as  the  exterior  walls 
but  was  composed  of  two  1/2-inch  sheets  and  a  3/4-inch  sheet. 


The  house  models  were  built  on  skids  enabling  them  to  be  towed  and 
placed  into  different  configurations  and  orientations.  In  addition,  single 
walls  of  some  of  the  houses  were  built  and  could  be  moved  from  place  to  place 
to  produce  additional  shielding  effects. 

5.2  FREE-IN-AIR  ENVIORNMENT. 

The  VCS  code  system  is  designed  to  couple  the  adjoint  calculation 
of  building  response  with  free-in-air  environments  calculated  with  the  dis¬ 
crete  ordinates  code,  DOT.  Calculations  with  DOT  of  the  free-field  radiation 

environments  produced  by  the  BREN  ^®Co  source  and  the  BREN  HPRR  reactor  used 
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in  this  study  have  been  performed  by  Dolatshahi  and  Kaul.  Here  we  provide  a 

summary  of  the  comparison  of  their  calculations  with  the  measurements  of  the 

FIA  kerma  made  during  the  BREN  experiments.  A  comparison  of  the  kerma  computed 

by  DOT  with  measurements  for  BREN  ®®Co  source  at  a  height  of  1125  feet  is  pre- 

2 

sented  in  Figure  2.  The  "r  kerma"  as  a  function  of  slant  range  is  presented 
at  the  air/ground  interface.  This  comparison  and  at  other  detector 
heights  shows  excellent  agreement  between  the  FIA  calculations  and  measure¬ 
ments.  Similar  comparisons  of  the  neutron  kerma  from  the  HPRR  reactor  are  pre¬ 
sented  in  Figure  3.  The  DOT  calculation  for  the  neutron  kerma  is  in  excellent 
agreement  with  measurements  at  all  heights  for  all  ground  ranges.  However, 
calculated  gamma-ray  kerma  associated  with  the  HPRR  source  do  not  agree  well 
with  experiment,  as  indicated  in  Figure  4  with  the  calculations  consistently 
low  by  about  40%.  As  discussed  in  Ref.  23,  the  cause  of  this  discrepancy  is 
unresolved,  although  it  is  possible  that  the  gamma-ray  detectors  were  not 
properly  calibrated  for  the  high  energy  gamma  rays  from  the  fission  products 
and  the  neutron- induced  capture  gamma  rays.  Neutron  interactions  with  the 
detector  system  itself  may  also  account  for  the  high  gamma-ray  measurement 
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XPOSURE  RT  (ROENTCENS/HR) 


SLANT  range  (K^") 


BREN  reactor  gamma-ray  kerma  rate  versus  slant  range 
measured  and  calculated  at  the  ground  surface. 


reported.  The  same  calculational  technique  has  been  applied  to  a  similar 
reactor  operating  at  the  Ballistic  Research  Laboratory  in  Maryland.  Compari¬ 
sons  of  calculated  garrma-ray  spectra  with  measurements  made  with  a  NE-213 


detector  show  much  better  agreement  than  the  integral  kerma  comparisons  for 
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BREN. 

All  our  house  transmission  calculations  use  the  energy  angular 
fluence  tape  produced  by  Dolatshahi  and  kaul's  DOT  calculations  of  the 
BREN  reactor  and  ^*^Co  experiments. 

5.3  calculations  FOR  HOUSE  A. 

The  first  set  of  house  transmission  calculations  were  performed  for 
a  house  model  identified  as  "House  A"  in  the  BREN  experiments.  House  A  was 
designed  to  replicate  a  standard,  one-story,  typical  Japanese  residential 
house  present  at  both  Hiroshima  and  Nagasaki.  A  computer-drawn  representation 
of  our  combinatorial  geometry  model  of  House  A  is  presented  in  Figure  5. 

The  House  A  geometry  model  for  the  radiation  transport,  as  well  as  other 
houses  used  in  the  calculations,  is  based  on  as-built  drawings  obtained  from 
Nevada  test  site  personnel  and  from  dimensional  measurements  made  in  1982  by 
the  authors  on  the  house  frames  that  were  still  standing  at  the  Nevada  test 
site.  Some  of  the  dimensions  for  House  A  are  indicated  in  the  figure. 

Samples  of  the  cement  asbestos  board  used  for  the  house  walls, 
partitions  and  roofs  were  obtained  from  the  Nevada  test  site  in  the  vicinity 
of  the  houses,  and  an  elemental  analysis  was  performed  at  the  Oak  Ridge 
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TYPICAL  DIMENSIONS 

Length:  28'  4” 

Width:  22'  4" 

Roof  Peak:  13'  6" 

Interior 

Openings 

Height:  5'  8" 

Interior 
Partitions 
Height:  9'  2" 


Figure  5,  House  A  geometry  model. 


/  / 
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National  Laboratory.  The  result  of  the  analysis  is  presented  in  Table  2. 

Calculations  for  the  ^*^Co  source  used  the  "PVC"  cross  section  set^^  consisting 

of  Pg  Legendre  expansion  of  the  scattering  angular  distributions  and  36  multi- 

groups  for  gatnna-ray  transport.  Calculations  for  the  HPRR  source  used  the 
26 

DLC-31  P^  expanded  cross  sections  in  coupled  multigroups  (37  neutron  group  - 
21  gamma-ray  group).  Elemental  cross  sections  available  in  both  libraries  (in¬ 
dicated  in  Table  2  )  were  mixed  at  the  measured  concentration  to  form  the  macro¬ 
scopic  cross  sections  for  dry  TRANSITE.  Those  elements  without  cross  sections 
compose  0.04%  of  the  mass  of  the  composition.  Our  calculations  employed  a 
partial  density  2.08  gm/cc  of  TRANSITE  plus  a  partial  density  of  0.023  gm/cc 

of  moisture.  Simple  geometry  ^*^Co  attenuation  measurements  of  the  house  ma- 
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terial  were  made  at  our  laboratory  fay  Shreve.  A  comparison  of  the  atten¬ 
uation  measurements,  with  the  attenuation  predicted  by  the  cross  sections 
used  in  our  calculations,  is  shown  in  Figure  6. 


Calculations  using  VCS  were  performed  for  House  A  for  several  dif¬ 
ferent  locations  and  detector  heights.  Our  calculations  for  House  A  and 
other  configurations  were  not  designed  to  be  exhaustive  but  a  representative 
comparison  with  the  variety  of  measurenents.  Furthermore,  we  have  not  edited 
our  results;  every  calculation  that  we  have  performed  to  compare  with  BREN 
experiments  is  presented.  The  detector  locations  are  indicated  in  Figure  7 
which  presents  a  computer  drawing  of  our  House  A  model  looking  down  from  the 
top  with  no  roof  and  showing  the  interior  partitions  and  openings.  The 
detector  identification  system  is  identical  to  the  one  used  to  report  the 
BREN  measurements.®  The  house  was  oriented  so  that  the  large  opening  (the 
veranda)  faced  the  source.  The  source  was  located  at  a  ground  range  of  750 
yards  and  height  of  1125  feet.  Comparisons  of  the  calculated  and  measured 
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transmission  factors  for  House  A,  for  the  °°Co  source  and  the  HPRR  source  are 
presented  in  Table  3.  For  both  measured  and  calculated,  the  transmission 
is  the  in-house  kernia  divided  by  the  corresponding  free  field  kerma.  The 
HPRR  gamma  includes  prompt  gamma  rays  emitted  from  the  reactor,  secondary 
gamma  rays  produced  in  the  air  and  ground  and  secondary  gamma  rays  produced 

in  the  house  materials.  The  comparison  between  the  measured  and  calculated 
data  for  the  ^*^Co  gatma-ray  source  is  excellent.  The  neutron  transmission  fac¬ 
tors  are  also  in  excellent  agreement  with  the  measured  data.  Since  part  of  the 
calculational  procedure  involves  the  numerical  calculation  of  the  free-in-air 
environment,  it  is  not  possible  to  provide  a  computed  estimate  of  the  statisti¬ 
cal  uncertainty  in  the  calculated  results.  However,  the  fractional  standard 
deviations  involved  in  the  adjoint  Monte  Carlo  calculation  are  all  less  than  5*. 

The  comparison  of  the  calculated  gamma-ray  transmission  from  the  HPRR 
source  with  the  measured  values,  however,  is  not  very  good.  At  the  two  interior 
locations  attempted,  the  calculated  values  are  consistently  lower  than  the 
measured  data.  The  contributions  to  the  gamma-ray  transmission  factor  calcu¬ 
lated  are  provided  in  Table  4.  All  values  in  the  table  are  normalized  to  the 
FIA  total  calculated  gamma-ray  kerma.  At  the  location  of  the  house  (750  yards 
ground  range),  approximately  75%  of  the  free-in-air  kerma  comes  from  secondary 

gatima  rays  produced  by  the  neutrons  in  the  air  and  in  the  ground,  and  the 

other  25%  are  from  the  prompt  gatmna  rays  emitted  from  the  reactor.  The  neutron 

kerma  at  this  point  is  almost  three  times  higher  than  the  ganina-ray  kerma.  At 

the  five-foot  detector  height  for  location  E,  our  calculation  of  the  gamna-ray 
kerma  is  nearly  identical  to  the  free-in-air  ganna-ray  kerma.  In  this  instance, 
the  gaitiTia  rays  produced  from  neutron  capture  in  the  TRANS ITE  material  offsets 
the  attenuation  of  the  incident  gamma-ray  kerma.  A  similar  situation  exists 
for  the  three-foot  location  H  detector. 


. .. ...........  .s .  ... . . . 


Table  4. 

Gamma  transmission  analysis  for  house 

A  and  APRR 

source . 

Location 

H  3FT 

Location  E  5  Ft 

Component 

FIA  Kerma 

Gamma  TF 

Gamma  Kerma 

Gamma  TF 

Gamma  Kerma 

Air/Ground 

Secondaries 

0.744* 

0.323 

0.240 

0.362 

0.269 

Prompt 

0.256 

0.441 

0.113 

0.488 

0.125 

Total  Gamma 

1.00 

0.353 

0.353 

0.394 

0.394 

Neutron 

2.90 

0.214^^ 

0.621 

0.215^ 

0.624 

TOTAL 

CALCULATED 

0.97 

1.02 

MEASURED 

1.32 

1.44 

All  values  normalized  to  the  calculated  FIA  total  gamma-ray  kerma. 
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Multiples  neutron  FIA  kerma  to  give  gamma  kerma  fraction  (e.g.,  0,621  at 
location  H,  3  ft)  from  secondary  gamma-ray  production  in  house  materials. 


Since  the  gamma-ray  transmission  factors  for  the  prompt  gamma  rays 
and  air  secondary  gamma  rays  appear  reasonable,  it  was  conjectured  that  the 
discrepancy  between  the  calculated  and  measured  values  stems  from  the  produc¬ 
tion  of  gamma  rays  in  the  house  materials.  To  test  this  supposition,  a  one- 
dimensional  spherical  discrete  ordinates  calculation  with  the  ANISN  code  was 
performed  for  a  shell  source  of  the  BREN  HPRR  neutron  and  ganna-ray  spectrum 
at  750  yards  ground  range  incident  on  spherical  shells  of  varying  TRANSITE 
thickness.  The  calculated  ganma-ray  transmission  factor  at  the  center  of  the 
sphere  is  given  in  Figure  8,  The  gamma  transmission  factor  is  a  ratio  of  the 
gamma-ray  kerma  at  the  center  of  the  sphere  to  the  gaimia-ray  kerma  at  the  shell 
source  position.  As  the  TRANSITE  increases  in  thickness,  the  ganma-ray  trans¬ 
mission  factor  increases  because  of  the  increased  production  of  capture  gamma 
rays.  However,  beyond  a  certain  thickness,  the  self-shielding  of  these  capture 
gamma  rays  in  the  TRANSITE  reduces  the  transmission  factor.  The  maximum  trans¬ 
mission  factor  that  we  could  calculate  by  this  method  was  1.2.  The  point  on 
Figure  8  is  plotted  at  the  wall  and  roof  thickness  of  TRANSITE  for  House  A  and 
indicates  the  transmission  calculated  by  Monte  Carlo  for  House  A.  The  trans¬ 
mission  is  less  than  the  spherical  model,  as  expected,  because  the  openings  and 
rectangular  geometry  reduce  the  total  capture  source  and  transport  to  the  detec¬ 
tor.  This  indicates  that  it  would  be  impossible  to  achieve  transmission  factors 
on  the  order  of  1.4  with  the  cross  sections  used  in  our  calculations  because 
this  spherical  model  maximizes  the  contribution  from  neutron  induced  gamma-rays. 
Since  excellent  agreement  is  obtained  for  the  neutrons  and  ®®Co  gamma-rays, 
the  discrepancy  does  not  lie  in  our  geometry  model  of  the  house  or  the  radia¬ 
tion  transport  procedure. 
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Figure  8.  Maximum  transmission  factor  model  calculation 
for  TRANSITE  material. 
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It  is  possible  that  a  strong  secondary  gamma-ray  producer  is  pres¬ 
ent  in  the  TRANSITE  so  that  either  the  cross  sections  used  in  the  computations 
for  the  given  composition  are  in  error,  the  large  producer  is  among  the  ele¬ 
ments  detected  but  not  incorporated  into  the  calculations,  or  the  large  pro¬ 
ducer  was  not  detected  in  the  analysis  of  the  TRANSITE  samples.  The  elements 
not  incorporated  in  the  cross-section  set  for  the  calculations  account  for 
0.04%  of  total  mass  (see  Table  2).  A  bounding  analysis  of  gamma  production 
for  these  elements  indicates  they  could  not  account  for  the  60%  increase  in 
the  neutron  to  gamma  transmission  factor  (see  Table  4)  required  to  produce 
agreement.  By  the  same  token,  it  is  possible  but  not  probable  that  the 
cause  is  due  to  cross-section  errors  or  missing  elements  in  the  analysis. 

Additional  evidence  indicates  that  the  ganma-ray  transmission 
measurements  may  be  in  error.  We  examined  similar  house  transmission 
factor  measurements  made  during  the  nuclear  weapon  explosions  of  Opera¬ 
tion  HARDTACK. During  operation  HARDTACK,  the  House  A  TRANSITE  houses 
were  exposed  to  low  yield  nuclear  weapons  with  a  hard  neutron  output  spectra 
similar  to  BREN  and  at  source  elevation  angles  similar  to  the  BREN  experi¬ 
ments.  A  comparison  of  results  reported  from  the  HARDTACK  measurements  with 
the  BREN  data  is  shown  in  Table  5.  We  also  present  our  calculations  for  the 
BREN  experiments.  Note  in  Table  5  that  the  neutron  transmission  factors  are 
reasonably  consistent  with  one  another,  but  that  our  calculations  of  the 
gamma  transmission  for  BREN  agree  better  with  the  HARDTACK  data  than  the 
measured  data  at  BREN.  The  ganma-ray  measurements  at  HARDTACK  were  made 

with  a  chemical  dosimeter  system.  In  conjunction  with  the  same  anomalous 

23 

high  gamma-ray  measurements  in  the  FIA  kerma,  we  conclude  at  this  time 
that  the  measured  gamma-ray  reported  doses  at  BREN  are  probably  in  error  due 
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Comparison  or  c  ransmiss  ion  data  for  house  A  at  3REN  and  H.ARDTAC’.i. 
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Neutron  Transmission _  Gamna  Transmission 


Location 

BREN 

Measured 

BREN 

Calculated 

HARDTACK 

Measured 

BREN 

Measured 

BREN 

Calculated 

hardtack 

Measured 

5'  E 

0.46 

0.44 

0.41 

1.44 

1.02 

1.08 

3'  H 

0.41 

0.43 

0.35 

1.32 

0.98 

0.98 

J 


to  either  a  calibration  error  or  to  neutron  contamination  of  the  Geiger 
detectors . 

It  is  instructive  to  examine  the  in-house  kernia  as  a  function  of 
the  point  of  incidence  of  the  external  radiation  field.  This  data  can  be 
invaluable  devising  generalized  transmission  factor  models  for  application 
to  Japanese  A-bomb  dosimetry.  Tables  6  and  7  present  the  fraction  of 
the  exposure  kerma  according  to  its  incidence  on  the  large  front  opening, 
front  and  rear  roof.  "Front"  and  "rear"  refer  to  directions  facing  toward 
and  away  from  the  source  respectively  and  do  not  correspond  in  this  case  to 
what  is  normally  called  the  "front"  and  "rear"  of  the  house.  For  example, 
in  Table  6,  0.2  x  (FIA  kerma)  was  detected  at  the  3  ft  B  location  from 
radiation  incident  on  the  house  through  the  front  opening.  The  residual 
column  tabulates  kerma  from  incidence  on  the  remaining  house  features.  Data 
in  Table  6  for  the  ®^Co  indicate  that  for  these  detector  locations  the  radi 
tion  through  the  front  opening  and  front  roof  dominates  the  detected  kerma. 
As  one  moves  further  back  into  the  house  the  reduction  in  the  transmission 
is  caused  by  the  reduction  in  the  kerma  from  radiation  incident  on  these 
house  features.  Similar  conclusions  can  be  made  from  the  data  exihibited 
in  Table  7  for  the  neutron  source.  The  higher  fraction  of  dose  from  the 
residual  components  indicate  the  neutron  environment  is  less  forward  peaked 
as  expected. 
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Table  7.  Fraction  of  t<.)tal  exposure  neutron  ke 
radiation  Incident  on  house  coraponent 
range  and  HPRR  source  at  height  of  11 


Of  this  residua]  0.10  is  due  to  radiation  incident  on  an  adjacent  window 


5.4 


CALCULATIONS  FOR  HOUSE  B. 


A  model  for  a  two  story  Japanese  house  used  at  BREN  was  called 
House  B.  Our  rendition  of  House  B  used  in  the  transport  calculations  is 
shown  in  Figures  9  and  10.  We  calculated  transmission  factors  for  the  HPRR 
neutron  and  ^^Co  sources  at  detector  locations  indicated  in  Figure  10.  The 
comparison  of  our  calculations  with  experimental  data  is  given  in  Table  8. 

The  source  ground  range  and  height  was  identical  to  those  for  House  A.  The 
source  location  column  in  Table  8  refers  to  the  direction  toward  source 
referenced  to  the  directions  indicated  in  Figure  10.  These  comparisons 
show  excellent  agreement  between  calculation  and  measurement.  The  geometry 
of  House  B  is  considerably  more  complex  than  House  A  and  the  good  agreement 
with  some  "difficult"  detector  locations  indicates  the  ability  of  the  trans¬ 
port  calculations  to  adequately  treat  the  effects  of  partitions,  windows  and 
multifloors . 

5.5  CALCULATIONS  FOR  MULTIPLE  HOUSE  ARRANGEMENTS. 

A  high  building  density  existed  in  typical  residential  neighbor¬ 
hoods  of  Hiroshima  and  Nagasaki  ATB.  In  these  built  up  areas  with  houses 
close  to  another,  the  shielding  effect  of  adjacent  houses  is  important.  To 
examine  the  multiple  house  shielding,  the  BREN  experimenters  placed  the 
houses  close  together  in  a  variety  of  arrangements  and  orientations  with 
respect  to  the  sources  of  radiation. 

Two  multiple  house  arrangements  were  calculated  in  our  validation; 
one  for  ^*^Co  source  and  one  for  the  HPRR  neutron  source.  The  BREN  experi¬ 
menters  did  not  use  an  identical  arrangement  and  orientation  of  House  A  and 
House  B  for  both  ^^Co  and  the  HPRR.  The  two-house  configuration  using  House 
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A  and  House  B  that  we  modeled  for  the  ^^Co  measurements  is  shown  in  Figure 
11.  In  this  configuration,  the  two-story  House  8  shields  the  single  story 
House  A  from  the  ^°Co  source  which  was  placed  at  the  standard  750  yards 
ground  range  from  the  houses  and  at  height  of  1125  feet. 

We  compared  calculations  with  experiment  at  three  detector  loca¬ 
tions  in  House  A.  The  results  are  given  in  Table  9.  The  agreement  is  very 
good.  Also  shown  in  Table  9  for  comparison  is  our  calculations  of  the  trans¬ 
mission  factor  for  House  A  without  the  intervening  House  B  (reproduced  from 
Table  3  ).  This  data  indicates  that  House  B  in  this  configuration  provides, 
roughly,  the  same  amount  of  shielding  as  House  A  containing  the  detector. 

The  configuration  used  for  the  HPRR  source  is  shown  in  Figure  12. 
Note  that  the  arrangement  of  House  A  with  respect  to  House  B  and  the  orien¬ 
tation  of  the  source  are  different  than  the  ^^Co  source  configuration. 

Three  locations  in  House  A  and  one  location  in  House  B  were  calculated  and 
the  results  are  provided  in  Table  X.  The  agreement  between  calculation  and 
experiment  is  excellent.  A  surprising  result  was  the  increase  in  transmission 
at  location  E  with  the  addition  of  House  B  on  the  side.  Evidently,  House  B 
acts  as  a  scatter  source  of  neutrons  which  increases  the  kerraa  in  House  A  at 
location  E.  The  situation  is  reversed  for  the  5'  H  location.  Note,  however, 
that  for  location  H,  House  B  blocks  a  large  opening  on  the  side  and  the  mid¬ 
room  partition  shields  scattered  neutrons  from  the  front  portions  of  the  side 
wall  of  House  B.  These  results  are  not  obvious  a  priori,  but  the  agreement 
between  experiment  and  calculation  for  both  the  single  and  multiple  config¬ 
uration  confirms  them. 


Table  9.  Comparison  of  calculations  with  measurements  of  transmission 
factors  for  the  ^^Co  source  in  a  multiple  house  configuration 


Detector 

Location 


Detector 

Height 


Multiple  House 
Confiouration 


Measured  Calculated 


House  A 
Standalone 
Calculated 


Table  10.  Comparison  of  calculations  and  measurements  ot  neutron 
transmission  factors  for  the  HPRR  source  in  a  multiple 
house  configuration. 


Detection 

Location 


Detector 

Heiaht 


Multiple  House 
Configuration 


Measured  Calculated 


House  A 
Stand  Alone 
Calculated 
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CONCLUSIONS 


The  BREN  house  transmission  experiments  provide  an  excellent  set  of 
measurements  to  validate  the  calculational  procedures  that  will  be  used  to 
derive  house  shielding  estimates  for  the  revised  dosimetry  of  the  survivors 
of  the  Hiroshima  and  Nagasaki  A-bombs.  The  BREN  experiments  were  performed 
in  realistic  full  scale  models  of  Japanese  residences.  Although  the  radia¬ 
tion  spectra  and  relative  intensities  of  neutrons  and  gamma  rays  incident 
on  the  houses  from  the  HPRR  and  the  ^^Co  sources  are  not  appropriate  for 
direct  application  to  the  A-bomb  survivors,  they  cover  the  full  energy  range 
of  importance.  The  codes  and  calculations  required  to  compare  with  BREN 
experiments  are  the  same  as  needed  for  the  A-bomb  dosimetry,  Thev  consist 
of  a  two-dimensional  discrete  ordinates  calculation  of  the  FIA  environment 
coupled  to  an  adjoint  Monte  Carlo  calculation  in  detailed  house  geometry. 

The  agreement  obtained  between  calculations  and  the  experiments  is 
excellent  for  neutrons  and  ®^Co  gamma  rays.  Every  house  transmission  calcula¬ 
tion  spanning  simple  to  complex  configurations  and  detector  locations  for  the 
^^Co  and  HPRR  neutron  cases  was  within  an  acceptable  margin  of  error.  The 
gariria-ray  transmission  calculations  for  the  reactor  source  did  not  agree  well 
with  the  experiments.  Analysis  of  this  discrepancy,  however,  strongly  indi¬ 
cates  that  the  problem  probably  does  not  reside  in  the  calculational  procedure 
but  in  the  measurements  themselves. 


In  conclusion,  we  believe  the  excellent  agreement  of  our  calcula¬ 
tions  with  the  BREN  experiments  validates  the  cal culational  procedure  we  plan 
to  apply  to  estimating  the  house  shielding  for  survivors  of  the  Hiroshima  and 
Nagasaki  A-bombs.  Certainly,  the  calculations  for  Hiroshima  and. Nagasaki  will 
involve  modifications  to  the  code  used  for  the  computations  reported  here, 
but  to  the  extent  that  these  modifications  involve  increased  calculational 
complexity  to  treat  more  realistic  materials  and  configurations,  the  bench¬ 
mark  established  by  these  comparisons  with  the  BREN  experiments  provides  a 
valuable  point  of  departure. 
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